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Abstract 

We study what might be called fractional vortices, vortex configurations with the mini- 
mum winding from the viewpoint of their topological stability, but which are characterized 
by various notable substructures in the transverse energy distribution. The fractional vor- 
tices occur in diverse Abelian or non-Abelian generalizations of the Higgs model. The global 
and local features characterizing these are studied, and we identify the two crucial ingredi- 
ents for their occurrence — the vacuum degeneracy leading to non-trivial vacuum moduli A4, 
and the BPS nature of the vortices. Fractional vortices are further classified into two kinds. 
The first type of such vortices appear when M has orbifold Z n singularities; the second 
type occurs in systems in which the vacuum moduli space A4 possesses either a deformed 
geometry or some singularity. These general features are illustrated with several concrete 
models. 



* Now at Theoretical Physics Laboratory at RIKEN since April 2009. 
t Now at the department of physics in Kyoto University since April 2009. 



1 Introduction 



Vortices appear in many different areas of physics, from fluid and plasma dynamics, solid-state 
physics, particle physics to cosmology. Usually certain topological properties lie behind their 
stability in time and in space. Typically, the energy distribution in the plane perpendicular to 
the vortex axis is peaked around the vortex axis, with a well-defined finite width in the vortex 
profile. This is certainly the case for the single (i.e. minimum- vorticity) type II Abrikosov- 
Nielsen-Olesen (ANO below) vortex [H [2], where the origin of the vortex stability lies in the 
first homotopy group, m(U(l)) = Z. As it turns out, however, when the gauge group and/or 
the matter content of the system are of more general kind than the standard Abelian-Higgs 
model (Landau-Ginzburg model) — U(l) gauge group and one charged scalar field — , a variety of 
interesting generalized vortex solutions appear. 

The present paper is concerned with a class of vortex-like solutions which might be called 
"fractional vortices". They are characterized by the minimal quantized vorticity (winding or 
magnetic flux) from the point of view of topological stability; nevertheless, their transverse 
profiles exhibit various non-trivial substructures as if they were made of smaller vortices, a little 
like a multi- vortex solution in the standard type II superconductor. But in contrast to the latter 
case the tension of each of the sub-peak is not quantized, and their relative weights, distances 
and shapes depend on the details of the system, such as the coupling constants, the scalar VEVs 
(vacuum expectation values) and the symmetry breaking pattern, etc. In all cases, a sub-peak 
cannot be removed by sending its position to infinity while keeping others in a finite region. 

The aim of this work is to show that such a phenomenon occurs very generally, and to make 
a preliminary study of these solutions, trying to find what characterizes the occurrence and 
substructures of these fractional vortices. As the arena of our study, we consider various Abelian 
and non- Abelian extensions of the Higgs model. The degrees of freedom will be a set of complex 
scalar fields with various charges and gauge fields. For definiteness and for simplicity, we take 
the models whose Lagrangians have the form of the bosonic sector of M = 2 supersymmetric 
gauge theories. They are natural generalizations of the Abelian-Higgs model in the Bogomol'nyi- 
Prasad-Sommerfield(BPS) limit. The constant which characterizes the VEVs and which forces 
the system into the Higgs phase, arises as the Fayet-Iliopoulos (FI) term in the supersymmetric 

: An analogous phenomenon occurs in a simple extension of the Abelian Higgs model — the Landau-Ginzburg 
model — with two coupled scalar fields. This can be realized in a certain unconventional superconductor [3]. On 
a broader prospect, our fractional vortices share also some features with the fractional instantons — the torons [4] 
or the calorons [5] which exist when the base space has a period such as a torus. Fractional vortices in a torus or 
a cylinder have also been studied [7] , see also [8] . 
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setting. Although most of our results are independent of supersymmetry, we shall mention also 
results more specific to the supersymmetric versions of the models, when appropriate. 

It will be shown that the fractional vortices can be further classified into two different classes. 
The first type exists when the vacuum moduli Ai has an orbifold Z n singularity. The second 
type occurs when the vacuum moduli A4 has a 2-cycle with a deformed geometry. 

This paper is organized as follows. In Sec. [2] we review the semi-local vortices in the extended 
Abelian-Higgs model (EAH). In Sec. [3] general properties of vortices in degenerate vacua are 
discussed. In Sec. H] we present various concrete models admitting fractional vortices based on 
a CP 1 target space. The first two models, in particular, provides the simplest examples of the 
fractional vortices of the first and second types, respectively. In Sec. [5] we discuss an SO(N) x U(l) 
gauge theory, which exhibits fractional vortices of the second type. Some useful technical details 
are collected in the Appendices. 



2 Semi-local vortex in the Extended Abelian-Higgs model 

Something quite non-trivial occurs already in the Abelian-Higgs model, if the number of charged 
fields is greater than one [TQl EEH EG] . The model is 

C = -^F ia ,F'» + V lt q(Trq?-±(qq1-tf , (2.1) 

where = <9 M — iA^ is the standard covariant derivative, q = (gi, g 2 , . . . , q^ t ) represents a set of 
complex scalar matter fields of the equal charge. This model is sometimes called the semi-local 
model since not all global symmetries, i.e. G = U (Nf) here, are gauged. Even if we restrict to the 
minimum vorticity, the vortex profile turns out to depend on the particular solution considered. 

In order to have a finite-energy (the energy per unit length — the tension) configuration, the 
complex scalar fields must asymptotically approach a vacuum configuration, 

M = { Qi }, ^| gi | 2 = e, (2.2) 

i=l 

far from the vortex center. By the SU(Np) global and U(l) local symmetry they can be chosen 
to be 

(q) = (g ,0,...,0), q = y/£, (2.3) 
breaking the global symmetry to SU(N{ — 1) x U(l) In other words, the vacuum moduli space is 

M = CP Ni ^ = SU(N t )/[SU(N{ - 1) x [/(!)] . (2.4) 
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The vacuum configurations M represent a non-trivial U(l) fibration over the vacuum moduli 
space CP f_1 . As the U(l) C U(Nf) part is gauged, its breaking does not lead to any further 
vacuum degeneracy. Since the first homotopy group of the vacuum configurations M ( 12. 2p is 
trivial, vortex solutions may not necessarily be stable. 

Indeed, for (3 = A/e 2 > 1 (i.e. type II superconductors), an ANO vortex solution embedded 
in the first flavor is found to be unstable against fluctuations of the extra fields (i = 2, 3, . . . , Nf) 
which increase its size; the vortex flux spreads out all over the transverse space [HI [12]. 

For (3 < 1 (i.e. type I superconductors), instead, an ANO vortex p],[2] embedded in one of the 
flavors is found to be stable. The origin of the stability of such a vortex can be traced to the fact 
that the asymptotic scalar field must be actually the vacuum configuration (12.31) modulo gauge 
transformations. The vortex winds a non-trivial U{1) fiber over the vacuum moduli CP f_1 : the 
relevant homotopy is 

7n(Z7(l)) = Z , (2.5) 

just as in the case of the ANO vortex (which indeed it is). 

In the interesting special (BPS) case, (3 = 1, we find a family of degenerate vortex solutions 
with the same tension, T = 2ir^. Except for the special point of the vortex moduli space (i.e. the 
space of solutions), which represents the ANO vortex (sometimes called a "local vortex"), the 
vortex has a power-like tail in the profile function, and the width of the vortex (thickness of the 
string) can be of an arbitrary sizcl. Far from the vortex center, the vortex configuration essentially 
reduces to the CP^ -1 sigma-model lump (or two-dimensional Skyrmion), characterized by 

7T 2 (CP^ 1 ) = Z . (2.6) 

In terms of an effective potential as a function of the vortex radius, the k = 1 (minimum- winding) 
sector of the system has the minimum at the origin for (3 < 1; at infinity for (3 > 1 (a "run-away 
vacuum" behavior); and has no potential — a flat direction — in the BPS case. 



3 Vortices in degenerate vacua 

There are two crucial ingredients which lead to the interesting varieties of degenerate vortex 
solutions in systems such as the extended Abelian-Higgs model (with (3 = 1) just considered: 

2 This type of vortex solutions has been termed "semi-local vortices" . Again although this is not an entirely 
adequate terminology we shall stick to it as it is commonly used in the literature. 
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Metric singularities 
Fig. 1: Vacuum moduli A4, fiber F over it, and possible singularities 

the vacuum degeneracy and the BPS saturated nature of the vortices. The first means that, in 
contrast to the cases in which the vacuum moduli is trivial (a point we must consider in general 
all vortex solutions in all possible points of the vacuum moduli simultaneously. See Fig. [TJ Even 
if we restrict ourselves to the minimally winding vortex solutions only — we shall do so in this 
article for definiteness — the vortices represent non-trivial fiber bundles over the vacuum moduli 
Ai. The BPS saturated nature of the vortices implies that the vortex equations linearize, and in 
turn half of the equations (the matter equations of motion) being solved by the moduli-matrix 
Ansatz, as is well known [20], EU [22] , see Eqs. ( 1A.5I) . ( 1A.6I) . The other equations (the gauge field 
equations) reduce, in the strong coupling limit or anyway sufficiently far from the vortex center, 
to the vacuum equations for the scalar fields. In other words, the vortex solutions tend to sigma 
model lumps. 

3 This is the case for the standard Abelian-Higgs model, of course, but so it is in the case of the U (N) gauge 
theory with N{ = N squarks in the fundamental representation. The latter model and its generalizations, after 
the discovery of the non- Abelian vortices in the color- flavor locked vacuum [13] [14] , has attracted considerable 
attention [15J-[36j. 
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3.1 Structures of the vacuum moduli 



Let us first consider what we regard as the global aspect of our vortices. More precisely, our first 
concern is the vacuum moduli Ai on each point of which the vortex solutions are defined. Let 
the symmetry group of the underlying system be 

K = L®G F , (3.1) 

where L is the local gauge group, while G F is the global symmetry group. Let M be the manifold 
of the minima of the scalar potential, the vacuum configuration M = {qi | q^T ! q = £ 7 }. The 
vacuum moduli At is given by the points 

P eM = M/F , (3.2) 

where the fiber F is the sum of the gauge orbits of a point in M 

feF = {q9 \ q9 = gq} ^ geL/L , (3.3) 

where we have taken into account the possibility that a given vacuum configuration might leave 
a subgroup Lq C L unbroken: 

L { q} = {l £L\£ q = q} . (3.4) 

In other words the vacuum moduli are made of the points of M in which gauge-equivalent points 
are identified. 

A subgroup of the global group 

G {q} C Gf , such that G {q} = {g f G G F | g f q = tq}, I G L , (3.5) 



represents the unbroken global symmetry group of the syste 




A vortex solution is defined on each point of A4, in the sense that the scalar configuration 
along a sufficiently large circle (S* 1 ) surrounding it traces a non-trivial orbit in F (hence a point 
in At). The existence of a vortex solution at a point / G F requires that 

7ri(F,/)^l; (3.6) 

a vortex corresponds to a non-trivial element of tti(F, f). The field configuration on a disk D 2 
encircled by S 1 traces Ai, apart from points at finite radius where it goes off M (hence from At). 

4 In other words, is the subgroup of Gf i.e. transformations which can be "undone" by — or equivalent 

to — a local gauge transformation. 
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Fig. 2: Far from the vortex center the fields trace, along a circle, the gauge orbits (/) regarded 
as a single point of M.. Inside the circle, the fields provide a map from R 2 (S 2 ) to A4. 



In other words it represents an element of ^(.M,]?) u, where p is the gauge orbit containing /, 
or 

P = *-(/) : (3-7) 

Ti is the projection of the fiber onto a point of the basis M.. The exact sequence of homotopy 
groups for the fiber bundle reads 

► 7T 2 (M, /) - 7T 2 (M, P ) - 7T! (F, /) -> 7^ (M, /) -> 7n (A4, J>) - - ■ • (3.8) 

where 7r 2 (M/F, /) ~ 7r 2 (.M,p). Note that in our application of such a sequence to the physical, 
vortex solutions, the reference point / or p appearing in the definition of the homotopy groups, 
corresponds to the field configurations along the large circle S 1 encircling the given solution, see 

Fig. m 

Given the points f,p and the space A4, the vortex solution is still not unique. Any exact sym- 
metry of the system (internal symmetry G^ as well as spacetime symmetries such as Poincare 
invariance) broken by an individual vortex solution gives rise to vortex zero modes (moduli), 
V. The vortex-center position moduli V ~ C, for instance, arise as a result of the breaking of 
the translation invariance in R 2 . The breaking of the internal symmetry (Eq. (13.51) ) by the 
individual vortex solution gives rise to orientational zero-modes in the U(N) models extensively 
studied in last several years. See [371 EH [39] for more recent results on this issue. 

5 This notion can be made more precise by considering the strong-coupling limit of our systems where the gauge 
field equations reduce to the vacuum condition, sec Eq. (|3.9p below. In the presence of an ANO like sub-peaks, 
however, this correspondence leads to a singular lump, as will be seen in several examples below. 
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Our main interest here, however, is the vortex moduli which arises from the non-trivial vacuum 
moduli M. itself. Due to the BPS nature of our vortices, the gauge field equation (see Eq. flA.3l) lPI 



Fl 2 = g] tfT'q - Z 1 ) , (3.9) 

reduces, in the strong-coupling limit (or in any case, sufficiently far from the vortex center), to 
the vacuum equation defining M. This means that a vortex configuration can be approximately 
seen as a non-linear a-model (NLcxM) lump with target space Ai (for non-trivial element of 
7r 2 (.M)). Various distinct maps 

S 2 i-> M , (3.10) 

of the same homotopy class correspond to physically inequivalent solutions; each of these corre- 
sponds to a vortex with the equal tension 



r, 



mm 



-e 7 J d 2 xF( 2 >0, (3.11) 



because of their BPS nature. They thus represent non-trivial vortex moduli. 

The semi-local vortices of the extended- Abelian Higgs (EAH) model reviewed in the previous 
section arise precisely this way. In the EAH model with N flavors of (scalar) electrons, 

M = S 2N ~ l , F = S 1 , M = S^/S = CP"- 1 , (3.12) 

and the homotopy sequence reads 

• • • - 7r a (S™- 1 ) - tt 2 (CP"' 1 ) - n t (S 1 ) - m (S 2 ^ 1 ) - . • . (3.13) 

1 Z Z 1 

The usual argument tells us then that ^(CP^ -1 ) and 7Ti(5' 1 ) are isomorphic: each (i.e. minimum) 
vortex solution corresponds to a minimal u-model lump solution. As in this model the vacuum 
moduli Ai is a (smooth) manifold, the above relations do not depend on the reference point / 
(or p). 

In most cases discussed below, however, the base space M. will be various kinds of singular 
manifolds: a manifold with singularities. The nature of the singularity depends on the system 
and on the particular point (s) of Ai. Some of them are simple conic (orbifold) singularities, due 
to the fact that some discrete (e.g. Z^) symmetry is restored at that point. The fiber (the gauge 



6 The index / denotes generally all the generators of the gauge group considered. A non-vanishing (FI) param- 
eter £ is assumed only for U(l) factor (s). 
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orbits) is smaller by some discrete quotient, with respect to F at neighbouring points. Other 
singularities at isolated points, or along some submanifold, reflect an even more drastic change of 
F such as a different unbroken gauge group at those points, as compared to that in surrounding 
regular (or less singular) points of Ai. The fiber itself goes through a discrete change in its 
dimension and in the type, at or along the singularity(ies). 



3.2 Classification of fractional vortices 



There are basically two distinct causes or mechanisms leading to the appearance of multiple peaks 
in the energy density even if the vortex under consideration has a minimum vorticity required by 
the regularity and the topological stability. The first type is related to the presence of orbifold 
singularities in Ai. For example, let us consider a TL 2 point po such as the one appearing in a 
simple U(l) model with two scalars (Section 14. ip . At this singularity, both elements of 712 (Ai,p) 
and 7Ti (F, f) make a discontinuous change. The minimum element of 7Ti (Fq, /q) is half of that of 
7Ti (F, f) defined off the singularity, and similarly for n 2 (Ai,po) with respect to 7i 2 (Ai,p), p ^ po- 
Even though the exact sequence such as Eq. (13. 8p continues to hold on and off the orbifold point, 
the vortex defined near such a point will look like a doubly-wound vortex, with two centers (if 
the vortex moduli parameters are chosen appropriately). 

Another cause for the appearance of fractional peaks, which we call the second type, is best 
understood by considering the strong coupling limit where the vortex reduces to a sigma-model 
lump, as already noted. Even if the base point p is a perfectly generic, regular point of Ai, not 
close to any singularity, the field configurations in the transverse plane (S* 2 ) trace a 2-cycle in 
the vacuum moduli space Ai. The energy distribution reflects the nontrivial structure of Ai as 
the volume of the target space is mapped into the transverse plane, C 

d 2 K 



E = 2 / xliZij W 18 ^ = 2 / BdK ■ ( 3 - 14 ) 

Jc d<j> l d<p J Jc 

Let us consider the case in which a 2-cycle in the vacuum moduli space Ai (the target space 
of the non-linear sigma model) is endowed with a deformed geometry, with regions of relatively 
larger scalar curvature, and possibly with some singularities. When the field configuration sweeps 
such regions, the energy density will show sub-peaks as illustrated in Fig. 



7 The existence of the directions in the target space, which are not related to any isometry, is necessary for the 
fractional lumps of the second type. Such directions are parametrized by so-called quasi-Nambu-Goldstone modes 
in the context of supersymmetric theories [40] while the directions of isometries correspond to Nambu-Goldstonc 
modes. 
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target space 
(vacuum moduli) 




deform 




lump configuration 




split 



Fig. 3: A sketch of fractional lumps of the second type in a sigma model limit. 



The field configuration may also simply hit one of the singularities (conic or not), which could 
represent a sick point in the non- linear sigma model limit. Even at finite coupling, the vortex 
tension density will exhibit a similar substructure. The existence of the singularity is, however, 
not essential for the occurrence of fractional vortices of the second type, in contrast to the first 
type. 



Several concrete models will be studied below. The fractional vortices appearing in these systems 
are caused by one or the other of the above mechanisms, or by collaboration of the two. The 
actual manifestation of these singularities could sometimes look quite complicated. We first 
discuss in this section models where the base space (vacuum moduli) is a CP 1 with one or two 
singularities, or a smooth but deformed CP 1 . 

4.1 Abelian Higgs model with two fields of different charges 
— fractional vortex of the first type — 

The first model is a simple extension of the Abelian-Higgs model with Nf = 2 flavors H = (A, B) 
but with unequal charges. We assign the U(l) charges ({m, n}) to the fields A and B, respectively. 
The gauge transformations take the form, 



4 Models based on CP 1 
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For simplicity, we assume that the charges are relatively prime, i.e., g.c.d({m, n}) = 1. The 
vacuum manifold (D-flatness condition) is topologically equivalent to S 3 and the vacuum moduli 
are topologically the same as CP 1 but with some conical singularities 

M = {A, B | m\A\ 2 + n\B\ 2 = £} , (4.2) 
M = M/U(l) ~ WCP^ n) ~ CP 1 /(Z m x Z n ) . (4.3) 

The vacuum moduli can be also described by the following quotient 

(A, B) ~ (X m A, \ n B) , A e C* . (4.4) 

Clearly, A = is a Z n fixed point and B = is a Z m fixed point. The t/(l) gauge symmetry is 
broken at every point of the vacuum moduli, thus topologically stable vortices can appear. 

Such a vortex solution is characterized by the broken ?7(l)-winding number v given in 
Eq. flA.4j) . The BPS energy density and mass are 



£ = + J=^HH^, (4.5) 

T = [ dx 2 S = 2tt& . (4.6) 



By using the moduli matrix method in Appendix, it can be expressed as 

v = -— [ dx 2 F 12 = - I dx 2 ddlog\s\ 2 , (4.7) 

2lT J 7T J 

W = -id logs, (4.8) 

H = (A, B) = (s~ m A (z) , s- n B (z)) , (4.9) 

where v is a positive number, s is an everywhere non-zero function and the moduli matrices A (z) 
and B (z) are polynomial functions of z. The first equation determines the asymptotic behavior 
of s as 

|s| 2 ^|^ as \z\^oo. (4.10) 

We choose the boundary condition 

(A, B) -> (A vcv e imue , B vcv e mue ) G M as \z\ oo . (4.11) 

The BPS equations (see Appendix) lead to the master equation 

e 2 

dd\ogu = -— [mu- m \A \ 2 + nuj- n \B \ 2 -Z] , (4.12) 

where to = ss^. 
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Before going into details, let us make a comment. If we fix A = (B = 0) everywhere, 
we can think of the system as just the Abelian-Higgs model with one complex scalar field B 
(A) whose U(l) charge is n (m). The vortices there are the normal ANO solutions, though the 
A;-vortex solutions will have the [/(l)-winding number k/n (k/m) with tension Tb = 2n^k/n 
{Tj\ = 2irl;k/m). Indeed, when only one field is active while the other is inert, the U(l) gauge 
coupling constant and the FI term can be rescaled such that the system looks exactly as the 
standard Abelian-Higgs model with unit U(l) charge. What we are trying to study in this 
section is an intermediate situation between two kinds of vortices where both fields contribute 
non-trivially. Such intermediate states should have the energy T = mT^ = nTs, and we shall 
see configurations which have m peaks in one limit and n peaks in another limit. 

First we choose a generic point such as A vcv ^ and B vev ^ 0. The moduli matrices behave 
asymptotically as follows 

A (z) = s m A -> \z\ mv e imve A wcw , B (z) = s n B -> \z\ nu e inw9 B vey , as \z\ -> oo . (4.13) 

Holomorphy of Ao,Bq requires mv G Z + and nv G Z + . As we have chosen m and n to be 
relatively prime, this is satisfied by v = k G Z + . Thus we have obtained the non-trivial condition 
for A ,B 

v = k : A {z) = A vcv z mk + 0(z mk - 1 ) , B (z) = B vcv z nk + 0(z nk ^) . (4.14) 

Note that k vortices have (m + n)k moduli parameters with the boundary vacuum modulus. 
They may correspond to positions and sizes of the fractional vortices. 

When we choose the special point A vcv = (Z m fixed point) or B vcv = (Z n fixed point) as a 



boundary condition, the conditions for the moduli matrix drastically change. Say \ A VCV \ = y^Jm 
and B vcv = 0. Immediately we get v = k/m and the conditions 

u= — : A = J^z k + ---, B = bz? + ---, (4.15) 
m V m 

where (5 is a semi-positive definite integer less than nv = —k. If we set Bq = 0, the solution 
is identical to the ANO vortex as we mentioned before. When Bq is not zero, the solutions 
significantly differ from the ANO solution and also from the semi-local vortices in EAH model. 



Similarly, if we choose \B YeY \ = \/£/n and A vcv = 0, the U(l) winding number becomes v = k/n 
and the conditions change as 



-: A Q = az a + ---, B = \^z k + ---, (4.16) 
n V n 
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where a is a semi-positive definite integer less than mv = —k. Note that the U(l) charge v is 
fractionally quantized at the conical singularities. The present model thus nicely illustrates the 
first mechanism for the fractional vortices discussed in the previous section. 

From this point on, we shall concentrate on the special concrete case m = 2 and n — 1, 
in order to illustrate in detail the properties of a fractional vortex. The vacuum moduli M. = 
WCPL jn ~ CP 1 /Z 2 has a Z 2 conical singularity at B vcv = (north pole), see Fig. HI We consider 



b = o 

\p = 00 




ip = 
A = 

Fig. 4: A sketch of the CP 1 /Z 2 which is a sphere with a conical singularity at a (i.e., north) pole. 

the minimal-energy vortex configuration (k = 1). When we choose a generic point (-B V ev 7^ 0) as 
the boundary condition, the minimal configuration has U(l) winding v = 1 whose energy is 

T min = 2tt£ , (u = 1) . (4.17) 

The corresponding moduli matrix is given by 

A (z) = A vcv z 2 + aiz + a 2 , -B (z) = B vcv z + b x , a u a 2 , h G C 3 , (4.18) 

with 2|A vev | 2 + |-B VCV | 2 = £. Although this is the minimal-energy configuration, we have three 
complex moduli parameters 01,02,61- Remember that A (B) is zero at a point where Aq (Bq) 
is zero. Note that Aq has two zeros and Bq has one zero because A winds twice and B winds 
once when we go around the boundary, S l at spatial infinity. An important observation is that 
the U(l) gauge symmetry is not generally recovered at the zeros. Only when A and B vanish 
simultaneously, the U(l) gauge symmetry is recovered (this would happen if some of the zeros 
of Aq and Bo are coincident). 
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Consider now the vortex at the special point of the vacuum moduli, B vev = 0. The minimal 
configuration k = 1 corresponds to v = 1/2 and has a tension 

TX cial = 7r£, (^ = 1/2). (4.19) 

The moduli matrix takes the form (for k — 1) 

A = A vcv z + a , 5 = 6 , A vev = , (4.20) 

where a, b are the moduli parameters. Comparing this with Eq. ( 14.181) with B vcv = 0, one 
immediately sees that the latter is not a minimal-energy solution. 

The vortex (energy, and magnetic flux) profiles can be approximately determined from the 
strong-coupling limit consideration. The gauge theory reduces to the non-linear sigma model 
whose target space is the vacuum moduli M. in Eq. (14.31) . The Kahler potential is given, in the 
supersymmetric version of our model, by 

K= \A\ 2 e~ 2V + \B\ 2 e~ v + £V . (4.21) 

Integrating out the U(l) vector multiplet V, we get the following Kahler potential in terms of 
an inhomogeneous coordinate: (f = 2y^A/B 2 

K = Zlogf(<p,<p)+£f- 1 (<p t <p) , /(p,£) = 1 + V1 + 2M 2 • (4-22) 

Note that the first term is due to the magnetic flux F\2 and the second term corresponds to the 
surface term d 2 J in Eq. (14.51) . All the regular BPS solutions are analytically solved by 

2C\s\ 2 = \B \ 2 + V\Bo\ 4 + 8£|A | 2 = |flb| 2 /fa <f>) , <P = <p(z) = ■ (4.23) 

Only the solutions which have points where A and B simultaneously vanish cannot be seen in this 
limit, because the U(l) gauge symmetry would remain unbroken there. Such solutions contain 
small lump singularities and we should go back to the original gauge theory in order to observe 
the configurations correctly. A numerical result is shown in Fig. 

As we move in the vacuum moduli space M. by varying the VEVs A vev , B vev (or ip vcv = 
2a/^t4 vcv / B 2 cv ) and change the vortex moduli parameters the tension density profile shows varying 
substructures. Since the zeros of the fields do not imply necessarily the restoration of a U(l) gauge 
symmetry, the positions of the peaks do not always coincide with the zeros of A, B. Although 
it is very complicated to specify the positions of peaks analytically, it is easy to visualize it 
numerically. In Fig. El we have shown the zeros of A, B and the peaks. We observe that there 
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Fig. 5: The energy (the left-most and the 2nd left panels) and the magnetic flux (the 2nd right panels) density 
arc shown, together with the boundary values (A, B) (the right-most panels) for the minimal lump of the first 
type in the strong gauge coupling limit. The moduli parameters are fixed as ai = 0, 0,2 = 1, 61 = —1 in Eq. (|4.18j) . 
The red dots are zeros of A and the black one is the zero of B. £ — 1. The last figures illustrates the minimum 
lump defined at exactly the orbifold point (see Eq. I|4.20[) ) with A vev = l/\/2, and with b = 0.8. 
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are no direct relations between the zeros of fields and the positions of the peaks, except at the 
two poles, A vcv = (south pole) and B vcv = (north pole), of the space M.. 

An axially symmetric peak appears at the zero z = z s of B (z) in the limit A vcv — > 0; as A vev 
departs from 0, it decomposes into two sub-peaks. We cannot remove one of the two sub-peaks 
pushing its position to infinity. This feature can be easily observed for large |y? VGV | 3> 1 where 
positions of the two peaks are naturally approximated by the zeros z = zf(i = 1,2) of Aq(z). 
(Here A vcv (zf + zf) = A vcv zfz2 = a 2 , B vcv z N = —bi.) The energy density E = 2ddK at 
those points is given by 



\zf 


- z s \ 


to 




-z N 
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E\ z=zf = 2fl»wr fj J.\, . (4.24) 



For instance, if the three zeros get separated by large distances, then we see that the sub- 
peaks are diluted. If, instead, only one of the zeros, z = is pushed toward infinity, that is 
\z% — z s \, \zf — z$\ ^> \zf — z s \, the peaks at z — zf becomes singular. In either case, the isolated 
one peak is not allowed as a vortex(lump) solution. This solution consisting the two sub-peaks 
is one of typical examples of fractional vortices. Only when B vcv = 0, they become independent. 
Such a limiting configuration is no longer a minimal energy configuration, however. The minimal 
configuration at exactly B vcv = (with (14.201) ) has only one peak. Its tension is half of the 
minimal configuration for B vcv ^ 0. 

The reason why the minimum vortex at B vcv ^ must have twice the energy with respect 
to the minimal object at B vcv = is as follows. Our vacuum moduli has a Z 2 singularity at 
|.B vev | = 0. If the vacuum is chosen at B vcv ^ the solution touches the singularity at a finite 
point in the z-plane and would get singular there. To remove such a singularity, the solution 
must wrap twice around the vacuum moduli. On the other hand, if one is at exactly the Z 2 
point the solution never touches it and a regular solution can be constructed with just a single 
winding. 

As discussed in Section [31 these characteristics of the vortex-energy profile are thus deeply 
rooted in the property of the vacuum moduli M. itself and to its singularity structure. In 
understanding the qualitative features of the vortex tension distribution and their dependence 
on p = ip vcv = 2^A vev / B^ cv just described, the crucial fact is that the elements of the homotopy 
groups corresponding to regular configurations make a discontinuous jump at p = ip = oo (-B V cv = 
0). In fact 

7T2 (M,p) Til (F, f) 

7r 2 (Al,oo) 7Ti(F,/ ) 
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The S 1 fiber itself reduces to half at the orbifold singularity 

/ = ^(p) = S 1 , /o = tT^oo) = S l /Z 2 . (4.26) 

Thus even though tt 2 (M) = 1 and n\(M) = 1, just as in the case of the EAH model Eq. (I3.13p . 
each minimum element of 7r 2 (A4,p) is a double cover of the minimum element of 7T 2 (A4, oo), just 
as the fiber at the generic p (S 1 : a = — > 2ir in Eq. (14.11) ) is a double cover of the fiber at B = 
(a = — > ir in Eq. (14.11) ). This is the (global) reason for the double peaks observed in Fig. [51 

The argument here can be easily extended to more general cases with the multiple flavors 
H = (A, B,C, D, ■ • • ) with generic U(l) charges Q = (m, n,o,q, ■ ■ ■), which are all relatively 
prime. The moduli manifold is then M = CP£~*, ~ CP Nf ~ 1 /(Z m x TL n x •■■). Near a 
Z m singular point, (|A vev |, |-B vev |, |C vev |, • • • ) = (y^Jrn, 0, 0, ■ ■ • , ) m peaks appear in the energy 
distribution. 

4.2 An Abelian U(l) x U(l) Higgs model 

— fractional vortex of the second type — 

The next system we consider, which has the same target space as in the previous model, is a 
U(l)i x U(l)2 gauge theory with three flavors of scalar electrons H = (A,B,C) with charges 
Qi = (2, 1, 1) for Z7 (l)i and Q 2 = (0, 1, —1) for Z7(l) 2 - The gauge transformations act as 

(A,B,C) -» (e^WAe^'^B^^^C) . (4.27) 

Note that the transformation (a, (3) = (ir, ±7r) leaves the fields invariant: the true gauge group 
is 

[U(l)x x U(l) 2 ]/Z 2 . 

Some details about this model are given in the Appendix. The vacuum manifold and vacuum 
moduli space are 

M = {A,B,C I 2\A\ 2 + \B\ 2 + \C\ 2 = ft, \B\ 2 - \C\ 2 = 6} , (4-28) 
M = Mj [(U(l)i x U(l) 2 ) /Z 2 ] . (4.29) 

Here ft is the FI term for the first U(l)i and ft is that for the second U(l) 2 . Note that the 
existence of the supersymmetric vacuum requires ft > |ft|. Below we shall mainly be interested 
in the case of ft = (in the Appendix the case with a non- vanishing ft is also discussed). 

When ft = 0, the vacuum manifold is the same as one in Section [4.11 In fact, the Kahler 
potential of the vacuum manifold is the same as Eq. (14.221) with replacements £ — > ft and 
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ip = 2y/2£A/ B 2 — > 2y/2£iA/(BC). Although the vacuum moduli manifolds are the same, there 
is an important difference. The singular point in the previous section was a Z 2 conical (orbifold) 
singularity whereas the singular point \B\ = \C\ = here represents a theory with a restored 
U(l)2 gauge symmetry, i.e., in a Coulomb phase. Since the Z 2 action has been modded out from 
the beginning, the singular point is not a Z 2 fixed point and therefore can be smeared by the 
introduction of non- vanishing £ 2 as illustrated in Figj6j (See however the discussions below.) 



Singularity (Coulomb phase) 




6 = & / 

Fig. 6: A sketch of the vacuum moduli space. 



BPS vortex solutions can be treated by the moduli matrix formalism. The tension is deter- 
mined by the f7(l)i winding number v\ (and the C7(l) 2 winding number v 2 if £2 7^ 0) as 

£ = -&F$ - + d 2 J , J = l -HH^ , (4.30) 

T = J dx 2 £ = 2tt^i + 27r£ 2 v 2 , (4.31) 
In the moduli-matrix formalism, the winding number can be expressed as 

vi = -^- [ dx 2 F$ = - f dx 2 dd\og \ Sl \ 2 , Wi = -id\og Sl , (4.32) 



27T J ±Z 7T . 

H = (A,B,C) = (si 2 A (z), s?s?B (z), s^s 2 C (z)) , (4.33) 

where s\ and s 2 are everywhere non-zero functions and the holomorphic functions Aq, Bo, Co are 
the elements of the moduli matrix. The BPS equation for the gauge fields reduces to the master 
equations 

oology = -j [wr 1 (2^r i |Aoi 2 + ^ 1 |5 | 2 + ^ 2 |c7 | 2 ) - a] , (4.34) 

dd\ogu 2 = [u^ 1 (w^lflol 2 - ^ 2 |C | 2 ) - 6] , (4-35) 



where u>i = Sis\, i = 1,2. 
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To specify a solution, we need to choose the winding numbers v\ and v 2 . The condition for 
v\ and v 2 depends on the boundary condition (A,B,C) — > (A vcv , B VCV ,C VCV ) e M as |z| — > oo. 
So we first need to determine (A vcv , B vcv , C vcv ) , (see the Appendix). It should be noted that, 

(2) 

when £2 = 0, although the magnetic flux F{ 2 does not contribute to the tension it contributes 

to the tension density through the surface term, J of Eq. (I4.30p . In other words v 2 as well as V\ 

(2) 

is needed to determine a solution. F 12 is non-trivial even when £ 2 = 0. 

Let us concentrate on the minimal energy configuration viz. £ 2 = in the following. We 
first choose a generic point (|B vev | = |C VCV | 7^ 0) as the boundary condition. The minimal 
configuration is given by two different choices of the winding numbers (the generic tension formula 
is given in Eq. ( 1A.19I) ) 

T min = <i , with {v u v 2 ) = (1/2, ±1/2) . (4.36) 

The corresponding moduli matrix is given by 

A (z) = A vcv z + a , B (z) = B vcv z + b , C (z) = C vcv , for (u u v 2 ) = (1/2, 1/2) , (4.37) 
A Q (z) = A vcv z + a , Bq(z) = B vcv , C (z) = C vcv z + c , for (v u v 2 ) = (1/2, -1/2) . 

The two complex parameters (a, b) represent the moduli parameters. 

At exactly B vcv = C vcv = the U(l) 2 gauge symmetry is restored at infinity: the system is 
in a Coulomb phase. We shall not discuss this case: it is beyond the reach of the moduli-matrix 
formalism. On the contrary, at an orbifold singularity such as those considered in the previous 
subsection the system remains in a Higgs phase even though with a different property from the 
surrounding vacua. 

Numerical solutions with various choices of A vcv , B vcv and C vev are shown in Fig. [71 When 
a point I-Bvevl = |C VCV | ~ is chosen as the boundary condition, one peak is observed near z a = 
—a/A vcv where A(z a ) = 0. On the other hand, another single peak appears near z b = —b/B vev 
where B(zf,) = 0, when we choose |A vev | = as the boundary condition. These two peaks 
are smeared in the intermediate values of (A vev , B vev , C vcv ). Although the vacuum manifold is 
the same as in Section 14.11 the vortex solutions are significantly different. The difference can 
be clearly seen in the winding numbers and tensions in Eqs. (14.1 7p and f!4.36p . In the model 
of Section 14.11 the Z 2 quotient requires us to choose v = 1 except for at the Z 2 fixed point. 
However, exactly in that point there are no such discrete symmetries, so there exist solutions 
with v\ = 1/2. On the contrary, in the model of this section, we cannot choose the singular point 
(in the Coulomb phase) as the boundary condition. 

In order to understand better the characteristics of the vortex energy profile in this model, 
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Fig. 7: The energy density (left-most) and the magnetic flux density (2nd from the left), (2nd from 
the right) and the boundary condition (right-most) are shown for the vortex of the second type. We have chosen 
& = 2, £2 = 0, ei = 1, e 2 = 2 and a = -1, 6 = 1 in Eq. (0T37]). 



we study the underlying sigma model, of which the Kahler potential with non-vanishing £2 is 
given by 



K = £if\(tp, ip) 1 + 61 log f x (tp, ip) - 1 2 2 log(l - A/ A (</?, v?)) , 
where /a(<£>, ^) is a function of an inhomogeneous coordinate tp = y/£iA/(BC) as 



1 - AM 2 + Vl + 2|^| 2 + A 2 |^| 



A = — . 



(4.38) 



(4.39) 



1+A ' " 6 

By setting £2 = (A = 0), we find the same Kahler potential as in Eq. (14.221) . Inside the disk of 
an arbitrary large radius the vortices reduce to sigma-model lumps, 

characterized by ir 2 (Ai). There the minimal lump solution has a tension Ti ump = 7r(£i — I^D- in 
the strong coupling limit giving rise to the above non-linear sigma model, solutions for s± and s 2 
are explicitly given by 



\BqCq 

V 1 - a/aO, <p) & 



S 2 



V 1 - W\(v,<p) 



5n 



Cn 



(4.40) 
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Since the vortex- lump solutions in the sigma- model limit, suffer from the Coulomb singularity, 
we introduce a small £2 to regularize such a singularity. For instance, the analytic solution 
I S2 1 2 = |A)|/|Co| in the sigma model limit (£2 = 0) leads to the singular magnetic flux for U(l)2 
since = —2dd log | S2 1 2 = —dd log \z\ 2 = —tt5^(z). The price we have to pay is that the two 
degenerate configurations (Eq. (I4.37P ) are now split. Suppose — £1 <C £2 < 0, then the energy 
changes as 

T (l/2,l/2) = 7r(6 + 6) = Tlump < T(i/ 2 -1/2) = 7t(£l ~ 6) , (4-41) 

although the lump solution cannot distinguish the two. The difference 27r|£ 2 | between the two 
is caused by the existence of an ANO-vortex-like singular peak attached on a tip of the lump 
solution for (z/i,z/ 2 ) = (1/2,-1/2). Actually, the solution for |s 2 | 2 with £ 2 < is still singular 
at the zero of Cq(z) since no configuration with C = can satisfy the D term condition for 
£2 < 0. Therefore, we consider only the vortex solution with (^1,^2) = (1/2, 1/2) in Eq. (14.371) 
as a smooth lump solution. The vortex configurations in a model with £2 in Fig. [8] (strong gauge 
coupling) are quite similar to the ones in the model with £ 2 = 0, Fig. [7| (finite gauge coupling). 
With (p vcv 7^ 0, 00, we again observe two peaks. The positions of the two are estimated by 
zeros of A (z) and B (z), z a = —a/A vcv , z b = —b/B vcv , which are mapped onto the two different 
poles of the vacuum moduli space Ai. Note that the two peaks for the fractional vortices in 
the previous subsection are mapped onto the same pole (the south pole ip vcv = 0). The energy 
density E = 2ddK at those zeros is 



-A 2 ) 


Vvev 


2 


2 


z 


a — Zb 


2 



1 - 


A| 


1 


A (^vev 


12 


z a — z b 


2 



n „ \r ,„ ^=6 ,;,, ' ,,_ _„ ■ (4-42) 



From this observation we see that the energy density for the fractional vortices gets diluted as 
\z a — Zb\ is increased. 

In spite of the similar structure of A4, the difference in the gauge group, matter content 
and in the fiber, manifest themselves in distinct ways that the vacuum manifold is covered by 
the vortex-sigma model lump solution here, as compared to the previous model. A comparison 
of Eqs. (14.171) and (I4.36P shows that the minimum lump solution in a generic point of Ai here 
clearly is seen to be half of the corresponding lump element ^(.M) in Section H~T1 (on the Z2 
singularity of the latter the solutions are similar). In the present model there is no jump in 
the homotopy-group elements at the singularity of M.. A sub-peak appears simply because the 
target space M. which is a distorted sphere is warped. If the target space were the standard CP 1 
sigma model (a perfect sphere) only one peak would have appeared. 

Just as the model considered in Section H~T1 showed a good example of the first mechanism for 



20 




Fig. 8: The energy density (left-most) and the magnetic flux density (2nd from the left), (2nd from 
the right) and the boundary condition (right-most) for the lump of the second type in the strong gauge coupling 
limit with a = -1 and b = 1 in Eq. (|P7|> . 6 = 2 and £ 2 = -0.1. 

the vortex substructures discussed in Section 13.21 the present model nicely illustrates the second 
mechanism for the fractional vortex. 

4.3 A model with U(l) x 577(2) 

The third and last example of a model with the target space of the droplet type is a gauge theory 
with a 17(1) x SU(2) gauge group with Higgs fields H = (A, B) 





C7(l) 


SU{2) 


A 


2 


1 


B 


1 


□ 
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namely, a complex scalar filed A of £7(1) charge 2, and two complex scalars B = (77i,7? 2 ) in 
the fundamental representation of the £77(2) group, and with the Abelian charge 1, while A is 
a singlet. The latter is conveniently denoted by a color-flavor mixed 2x2 matrix, B. There is 
furthermore a global symmetry £77(2)f. The gauge group acts on the fields as 

(A, B) -> {e 2ia A, e ia gB) , e ia e U(l) , g e £77(2) . (4.43) 

Note that the transformation (a,g) = (±7r, — 12) is a symmetry, thus the gauge group is really 
[17(1) x £77(2)]/Z 2 ~ U{2). 

The vacuum manifold is given by the JJ-flatness condition 

M = {(A,B) I 2|/4| 2 l 2 + 2BB t = £1 2 } . (4.44) 

By using SU(2) gauge and flavor symmetries, the vacuum configuration for B fields can be taken 
in the form B = £>1 2 , showing that the vacuum manifold M is an £ 3 defined by |^4| 2 + \B\ 2 = £/2. 
The vacuum moduli space A4 is a U(l) quotient of this £ 3 with the weighted charges given above, 
topologically the same as CP 1 . At a generic point (B 7^ 0) of the vacuum moduli, the flavor 
symmetry £C/(2)f is broken but the color-flavor diagonal symmetry £77(2) c +f is preserved. At 
the special point B = 0, neither the £77(2) gauge nor the flavor symmetry is broken. To be 
precise, we can compute the Kahler metric on the vacuum moduli space by eliminating all the 
gauge multiplets from the Kahler potential 

K = \A\ 2 e~ 2V + Tr ^BB t e -y-y '^ + £V , (4.45) 

with V being a U(l) vector multiplet and V being an SU(2) vector multiplet. Once we eliminate 
V from the Lagrangian, we obtain 
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K= \A\ 2 e- 2V + (2detB)2 e ~ v + £V . (4.46) 

Comparing this with Eq. f)4.2ip . it is easy to see that the two spaces have the same metric, if we 
replace the inhomogeneous coordinate tp = 2A/B 2 — > Aj detB. Thus the vacuum moduli space 
is of the droplet form as in Fig. HJ but there is no conical singularity. Instead, the SU (2) gauge 
symmetry is unbroken at the tip of the droplet (B = 0). 

Construction of the BPS vortex solutions in this model is summarized in the Appendix in 
detail. The energy formulae are the same as Eqs. (14.51) and (14.61) with J = \A\ 2 /2 + Tr (BB^)/2 
and the solutions are determined by U(l) winding number v 

v = -— [ dx 2 F 12 = - [ dx 2 dd\og\s\ 2 , (4.47) 

2ix J IT J 

W U{1) = -id log s , W su{2) = -iS'~ l dS' , (4.48) 
H = (A, B) = (V 2 Ao(^), ^ 1 £ / ~ 1 B (^)) , (4.49) 
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where s G C* and S' G SL(2, C) and Aq and all the elements of 'Bq(z) must be holomorphic 
functions of z. The BPS equations for the gauge the gauge fields lead to the master equations 

ddlogu = -j H-Tr (BoBjfi^ 1 ) + A| Ao |2 _ ^ 
8 (fiW" 1 ) = £ ^BoB^'- 1 - ^Tr (BoB^'- 1 )^) , (4.50) 

where uj = ss^ and Q! = S'S'\ 

The first equation determines the asymptotic behavior of \s\ 2 ~ \z\ 2u as \z\ — > oo. The 
choice of v must be consistent with the given boundary condition (A, B) — > (A vev , i3 vcv l2) satis- 
fying |A vev | 2 + |£> vcv | 2 = £/2. To find consistent solutions, it is useful to consider the following 
SU{2) gauge invariant / = detB = s _2 detB (2). When a generic point (B vcv 7^ 0) is chosen, 
holomorphy of Aq(z) and detBo(^) requires 1v G Z + because of the asymptotic behavior 

A ~ A vcv z 2 " + ■■■ , det B ~ B 2 ^z 2v + • • • , as |«| -»■ 00 . (4.51) 

In this way the vortices in this system are characterized by the half quantized U(l) winding 
number v G Z + /2. 

The minimal configuration with tension T = 7r£ is described by the moduli matrices 

A = A vcv z + a , det B = B 2 cv (z + b) . (4.52) 

Note that the matrix Bo (2) is not uniquely determined by these conditions. The simplest one is 
Bo = i3 vev diag(2; + 6, 1). This matrix breaks the color-flavor symmetry SU(2) c+ { of the vacuum 
into {7(l) c+ f. Henceforth, the generic configurations are generated by SU(2) c+ {, so that the 
vortex has an internal orientation CP 1 ~ SU{2) c+l /U{\) c+ { (Nambu-Goldstone mode). The 
moduli space of the single vortex is 

V = C x v2? , Vk3 = C x CP 1 , (4.53) 

which is a product of the orientational zero-modes (b) and a center of mass and a "size" parameter. 

When we choose A vev = as the boundary condition, the vortex is a semi-local extension of 
the so-called U(2) local vortex. In fact, if we turn off the size moduli, viz. a = 0, it is precisely 
the U(2) local vortex. The transverse width of the vortex becomes large as \a\ increases while 
the orientational moduli are still localized near the vortex core. This is completely different from 
the so-called non-Abelian semi-local vortex in U(N C ) gauge theories with iVf > N c Higgs fields 
in the fundamental representation. 
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We see that the global features of the vortices of this model are similar to that of the U(l) x 
£7(1) model (except, of course, for the internal, orientational modes which are present here). The 
key fact is that the non-Abelian sector of this model (with the matter fields B) is actually a 
U(2) ~ SU(2) x f/(l)/Z 2 theory. Taking into account the U(l) charge of A, we see that the fiber 
is generated by 7r 1 (S' 1 /Z 2 ), on and off the orbifold singularities (B = 0). 

The system at £> vcv = is in a (non-Abelian) Coulomb phase classically and no vortex 
solutions exist. In the U(l) 2 model of Section I4T21 we avoided the problem by turning on the FI 
parameter for the second U(l) factor. The same cannot be done here, as no Fl-like term exists 
for a non-Abelian gauge group. Therefore there exist neither a smooth non-linear sigma-model 
limit nor any regular lump solutions. Such solutions would necessarily be singular. 

4.4 An alternative U(l) 2 model: the lemon space 

The model of Section 14.21 was chosen to have a minimal fiber F with the fewest possible fields. 
However, the energy density (in the lump limit) and the Kahler metric turned out to be rather 
elaborate. In this subsection, we analyze a related model, with a slightly different charge assign- 
ment, Qi = (1, 1, 1) for Z7(l)i and Q 2 = (0, 1, -1) for U(l) 2 as 

(A, B, C) -> (e ia{x) A, e i"<P)+ift*)B, e" (l) -^ (l) C) . (4.54) 

An important difference from Section [4.21 is that the gauge symmetry is now t/(l)i x U(l)2, 
without a Z 2 division. The points (A,B,C) and (— A, B, C) related by (a, (3) = (tt, ±7r) G Z 2 , 
are distinct points. The vacuum manifold and vacuum moduli space are given by 

M = {A,B,C | \A\ 2 + \B\ 2 + |Cf =£i, \B\ 2 - \C\ 2 = , (4.55) 
M = M/ {U{l)x x U{\) 2 ) ■ (4.56) 

We see that A = is a Z 2 orbifold point, whereas the point B = C = represents a system 
in Coulomb phase (which can be Higgsed and regularized by £ 2 ^ 0). See Fig. [9j Clearly this 
model shares aspects both of the simple U(l) model of Section I4TT1 and of the U(l) x U(l) model 
of Section 14.21 In the following we shall consider mainly the case of £ 2 = 0, except when we 
consider the sigma-model limit, which is well defined only for a non-vanishing £ 2 . 

The vortex Ansatz is 

(A, B, C) = (s^A (z), s rVA)(*), s?s 2 C (z)) , (4.57) 
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Singularity (Coulomb phase) 




Z 2 singularity Z 2 singularity 

6=0 6^0 



Fig. 9: The lemon space. 

with the gauge field equations 

dd\ogu x = -j [lu^ 1 ({Aof + uj^lBof + LU^Col 2 )-^] , (4-58) 
^log^ 2 = --^^r 1 (^ 2 - 1 | J B | 2 -^ 2 |C , | 2 )-6] , (4.59) 

where u>i = s^sj, for i — 1,2. In order to avoid repetition, all the details are summarized in the 
Appendix. As in the case of Section I4~2"l the winding numbers are v\ and v 2 for C7 (l)i and £7(1)2, 
respectively. The tension depends only on v\ for £ 2 = 0: 

T = 2tt£iZ/i , i^i 6 Z + . (4.60) 

The minimal-energy solutions with the generic boundary condition (0 < |A VCV | 2 < £1) have 
T = 27r£i and are obtained by the following three different moduli matrices 

A = A vcv z + a, B = B VCV , C = C vcv z 2 + c x z + c 2 , (u ± , v 2 ) = (1, -1) , (4.61) 
A = A vcv 2 + a , .Bq = fi vcv 2; + b , C = C vcv z + c , {y u v 2 ) = (1, 0) , (4.62) 

A = A vcv z + a , B = B vcv z 2 + biz + b 2 , C = C vev , (ux, v 2 ) = (1, 1) . (4.63) 

As they obey different boundary condition for z/ 2 , they belong to different topological sectors. 
Each configuration has three moduli parameters. 

Near the Z 2 orbifold point we observe two peaks. Although the energy density always looks 
the same, the magnetic fluxes, especially of the second U(l) 2 , depends on the value of v 2 . In 
Fig. [TOl we show several numerical solutions for Eq. (14.631) . We also show a couple of solutions 
for Eq. (I4.62p in Fig. [TTJ In almost all regions, the configuration consists of one peak or two 
peaks but sometimes we observe three peaks simultaneously. 
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Singularity (Coulomb phase) 
(\A vev \\B vev \\\C ve 




-- (0.9,0.05,0.05) 



Z-2 singularity 

Singularity (Coulomb phase) 




(|A vev | 2 ,_B vev | 2 , |C vev | 2 ) - (0.7,0.15,0.15) 



I>2 singularity 

Singularity (Coulomb phase) 




= (0.4,0.3,0.3) 



1,2 singularity 

Singularity (Coulomb phase) 




j-BwIMCvbvI ) = (0.2,0.4,0.4) 



1>2 singularity 

Singularity (Coulomb phase) 




<{\A vev \\B vev \\ |C vev | 2 ) = (0.01,0.495,0.495) 
1>2 singularity 



\. 10: The energy density (left-most) and the magnetic flux density (2nd from the left), (2nd from 
right) and the boundary condition (right-most) for Eq. (|4.63[) with £i = 1 and £2 = and ei = 1, ei = 2. 
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Singularity (Coulomb phase) 




Z 2 singularity 

Singularity (Coulomb phase) 




1.2 singularity 

Fig. 11: The energy density (left-most) and the magnetic flux density (2nd from the left), (2nd from 
the right) and the boundary condition (right-most) for Eq. (14.62[) with £1 = 1 and £2 = and e± = 1, e<i = 2. 

On the other hand, at exactly a singular vacuum A vev = (the singular point on M), the 
minimal vortex with tension T — 7r£i is given by 

A = a, B = B vcv z + b, C = C VCV , {v x , v 2 ) = (1/2, 1/2) , (4.64) 
A = a, B = B VCV , C = C vcv z + c, {u u u 2 ) = (1/2, -1/2) . (4.65) 

At A = (</> = 0) a Z2 symmetry remains unbroken which is a typical orbifold singularity. As 
a result, the U(l)i fiber F is the half (a = — > tx) at the orbifold point as compared to that in 
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other points of the vacuum moduli, where a = — ► 2ir. The global structure of the vortex-sigma 
model lumps in this model is thus somewhat similar to the model of Section 14.11 At the Z 2 
orbifold singularity 7Ti(F) and 7r 2 (.M) make a jump, and this explains the appearance of the 
double peaks. 



As in the model in Section fl~2l we cannot take B vcv = C vcv = as a boundary condition since 
the second U(l)2 is unbroken at infinity. 

These aspects can be made more explicit in the strong gauge coupling limit ei,e2 — ► oo, 
where the Kahler potential is simpler (than in the model of Section 14.21) by construction and the 
solutions can be analytically solved. Since the Coulomb phase leads to singular solutions, we 
here turn on the another FI parameter £2 (IC2I < £1) f° r U(l)2- Working in a supersymmetric 
context, elimination of the gauge superfields V\ and V2 from 

K = \A\ 2 e~ Vl + \B\ 2 e- Vl - V2 + \C\ 2 e~ Vl+V2 + ^ + &V 2 , (4.66) 

yields the Kahler potential 

K = h log (l + V / A^+FW) " 161 log (|A| + V / ^TRS) . (4.67) 

where the inhomogeneous coordinate (p and A are defined by (p = 2 ^- and A = ^. The BPS 
solutions are given by the holomorphic functions 

m = , (4.68) 

and characterized by the quantized tension 

T = 2fV^, Ui = - ( dx 2 dB\og\ Sl \ 2 , (4.69) 

where |sj| 2 is given by 



= ( X + V^+O^W) , (4-70) 

2 \A \ 2 -A + y/A^ + (1 - Agg _ |j? | 2 2(1 - A) 

121 IC0I 2 2(1 + A) " |A | 2 A+ v/A 2 + (l-A 2 )|y5| 2 ' l ' ' 

4.5 An alternative U(l) x 5?7(2) model 

The next example with the same base space is the U(l) x S'C/ (2) theory with the same Lagrangian 



as in Section |4T3| except for a different U(l) charge assignment for the A field: 
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(7(1) 


dU[2) 


A 


1 


1 


B 


1 


□ 



The gauge group action on the fields is 

(A, B) -> (e ia 7L, e iQ #B) , e ia e U(l) , g G SC/(2) . (4.72) 

Note that (a,g) = (±7r, — 1 2 ) is not an identity operator, so that the gauge group is truly 
£7(1) x SU(2). The vacuum manifold is given by 

M = {(A,B) | |A| 2 1 2 + 2BB f = £1 2 } . (4.73) 

By using £77(2), we can bring B = £>1 2 . Then we see that the vacuum manifold is isomorphic 
to S 3 and the vacuum moduli space is isomorphic to S 3 /U(l) ~ CP 1 . The vacuum manifold is 
the lemon space as drawn in Fig. [9] There are two singularities: one (at A = ) is a Z 2 conical 
singularity and the other is a Coulomb singularity where SU{2) gauge symmetry is restored 
(B = 0). 

The construction of the BPS vortex in this model is the same as the one in Section 14.31 
However a difference appears in Eq. (I4.49P as 

H = [A, B) = (s^Aoiz), s-'S'-'Bo^ . (4.74) 

This makes a difference for the choice of the U{1) winding number v which should be chosen to be 
consistent with a given boundary condition (A, B) — > (A vcv , £> vev l 2 ) satisfying |A VGV | 2 + 2|£> V0V | 2 = 
£. To see this, we again make use of the holomorphic SU(2) gauge invariant / = det B = 
s~ 2 det Bo (z). For a generic point A vcv ^ or £> vcv ^ 0, the asymptotic behavior is of the form 

A -> A vcv z» + ■ • • , det B ^ Bl cv z 2v + ■■■ . (4.75) 

Holomorphy requires v to be semi-positive integer. Thus the minimal configuration with v — 1 
has the mass T = 27r£ and it is generated by the moduli matrix 

A Q = A vcv z + a , det B = Bl cv (z 2 + b x z + b 2 ) . (4.76) 

These conditions do not uniquely determine the matrix B (-2). The moduli matrix is the same 
as one for U(2) local vortex, for instance 

B (z) = B vcv ( 1 2 Cl [ +C \ ) • (4.77) 
z 2 + biz + b 2 
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Thus the moduli space for the single vortex is 

V gcn = C x V U k S , (4-78) 

where stands for the moduli space of k = 2 U(2) local vortices and the first factor C 

corresponds to the complex parameter a in A . Note that the minimal configuration includes 
two non-Abelian vortices since the vacuum moduli space has a Z2 singularity. Once we choose a 
generic point as the boundary condition, we cannot remove one of two non-Abelian vortices from 
the configuration. If we do that, the configuration meets a singularity. 

Only when we choose the special point A vcv = 0, we can avoid the Z 2 conical singularity. 
The U(l) winding number v can be a half integer and the minimal configuration is obtained by 
v= 1/2 

A = a , det B = B 2 vcv (z + b) . (4.79) 

This reflects the fact that at the orbifold point A = of A4, the U(l) fiber makes a jump 
(becomes a half of what is at regular points). The vortex moduli space is the same as that of 
k — 1 U(2) local vortex [T3J [TJ] and is given by 

V sp = Vfi? , v£? = C x CP 1 . (4.80) 

As in Section 14.31 we cannot take the singular point B vcv = as a boundary condition for 
constructing vortex solutions. The lumps in the strong gauge coupling limit always hit the 
singularity. 



5 U(l) x SO{N) model 

We now consider the fractional vortices occurring in a model with gauge group U(l) x SO(N). 
Vortices with orientational modes (non-Abelian vortices) in these models, in a maximally color- 
flavor locked vacuum, have recently been constructed and studied (37J [3H [39] . 

For our purposes here, we shall consider only the even-dimensional orthogonal groups, i.e. N = 
2M. The matter content is iVf = N flavors of squarks in the fundamental (vector) representation 
of the SO(N) group, all with the same unit charge with respect to the U(l) group: 





U{1) SO(N) 


H 


1 □ 
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As the Z 2 element (i.e. —1) of the SO(N) group is also an element of U(l), the gauge group is 
really 17(1) x SO(N)/Z 2 . 

The vacuum moduli have been studied in Ref. [JT] and it turns out that it has a rather rich 
structure. By color and flavor transformation, the scalar VEV can be put in the canonical form 

2M 

(H) = diag (vi, v 2 , •■• , v 2M ) , ^ v] = £ , ^el. (5.1) 

i=i 

Note that, in contrast to the U(N) models with Nf = N flavors, where vacuum conditions force 
the VEV of H to be proportional to an N x N unit matrix, the weaker condition here leaves the 
possibility of having arbitrary values Vi subject to the constraint, Y^h=i v f = £■ A l ar g e vacuum 
degeneracy is present here. 

At a generic point in At, where Vi ^ 0, Vi, and all distinct, the gauge and flavor groups 

L ^ uwxsow) Gr = su{2M)i (5 . 2) 

are completely broken. The fiber F is given by the L orbits of the points H = diag(t> i, v 2 , . . . , v 2 m) ■ 

On the points where some (at least two) of the v^s vanish, the unbroken gauge group Lo is 
strictly smaller than L. The gauge orbit F is now generated by L/L and has a smaller dimension 
than in the case of a vortex constructed on a generic point of M.. Thus even though in all cases 

7Ti(F) = Z , (5.3) 

its actual (e.g.) minimal element goes through discontinuous changes whenever we hit a singu- 
larity (or a singularity curve) on Ai. Also, in such a point, the global symmetry group Gp is 
different from that at surrounding points, and the consequent internal vortex moduli also under- 
goes a discontinuous change. As a singular surface (e.g. with a given number of vanishing v^s) 
contains a smaller subspace of singular points (some of the remaining v^s vanishing there), etc., 
one ends up with a rather rich structure of a (stratified) singular manifold Ai, and of the vortices 
and related sigma- model lumps as the fiber defined over it. We shall leave the study of these 
varieties of phenomena for a separate investigation: here we will take all vacuum moduli to be 
non-vanishing. 

Our fractional vortex solution is closely related to the "fractional lump" which was found by 
some of us recently [JT] • We choose in the following the scalar VEV to be of color-flavor diagonal 
form, and moreover proportional to the unit matrix form, 

(H) = v/ei 2 M (5.4) 
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leaving a residual global color-flavor symmetry SO (N) c+ { unbroken. The standard moduli-matrix 
Ansatz is 



H = s~\z, z)S'~\z,z)H {z) 



where s G t/(l) L , S' G SO(N) L . The gauge field BPS equations lead to 

.2 / 1 



dd log oj 

d ( n'dn'- 1 



AN \ oj 



Tr ( n n 



1-1 



where uo = ss^ , Q' = S'S'\ Qq = HqHq. J is the invariant tensor of SO(N = 2M) 



(5.5) 

(5.6) 
(5.7) 



J 



1m 
1m 



(5. 



The tension of the vortex remains 



T = 2v 2 / d 2 x ddlogu = 7rv 2 k . (5.9) 
./: 

Following the construction of [38], we have the constraint 

H^JH = z k J + O (z k - x ) , (5.10) 
for vortex solution of winding number k. 

In order to study the minimal winding vortex configuration more concretely, we choose 
zl M - Z C 



Z = diag(z!, z 2 , ■ ■ ■ , z M ) , C = diag(ci, c 2 , . . . , c M ) ■ (5.11) 



1 M/ 

To solve the master equations (15. 6p and (15.71) . we set 

Q' = diag ^e*' 1 , . . . , e x ' M , e - *' 1 , . . . , e" x « 
where the determinant one is manifest. Taking uo = e*, we obtain 



ddip 



8M 

2 



' M 

£{(l*-*l a + l«l a ) 

t=l 



ddx', = -- I (|* - Zi? + M 2 ) e-«**n - e-W-rt 



Vz G [1,M] 



(5.12) 



(5.13) 



(5.14) 



If we now take the infinite gauge coupling limit e — > oo, g — ► oo, we obtain the following lump 
solution 



e Xi = — Zi| 2 + |c- 12 



M 



M 



i=l i=l 



2 



(5.15) 
(5.16) 
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which has the energy density 

£ = 2£dd log | J] sf\z - ztf + |c,| 2 | . (5.17) 

This is the fractional lump solution found in Ref. [4T] . 

The vortex energy profile in the strong-coupling approximation for the U(l) x 5*0(6) model 
is shown in Figure [T2l Three fractional peaks are clearly seen. The positions of the peaks can 




-2 2 O.o' -2 2 



Fig. 12: The energy density of three fractional vortices (lumps) in the U(l) X 50(6) model in the strong coupling 
approximation. The positions are zi = —V2 + iv2, zi = — V% — Z3 = 2. Le/t panel: the size parameters are 
chosen as ci = C2 = C3 = 1/2. Right panel: the size parameters are chosen as ci = 0,C2 = 0.1, C3 = 0.3. Notice 
that one peak is singular (zi) and the other two are regularized by the finite (non-zero) parameters 02,3. 



be understood as follows. If q = one of the t/(l) C U(l) x SO(2M), constructed as the 
diagonal combination of U(l) and one of the U(l) Cartan subalgebra of SO(2M), is restored at 
the points z = Zi (i — 1,2, . . . M). The sharp peak in the right panel of Fig. [12] can be thought 
locally (in z) to be an ANO vortex. If q 7^ the situation around a fractional peak at z = z^ is 
more similar to the power-behaved semi-local vortex of the EAH model. The number of peaks 
reflects obviously the rank of the group considered (here rank{S'0(6)} = 3), but the number of 
the possible fractional peaks depends on the point of the vacuum moduli (a particular VEV) 
considered. For instance, if two of Vi are taken to be zero, the maximum number of the fractional 
peaks would be two, and so on. 

In the supersymmetric version of the models based on the U(l) x SO(N) gauge groups, the 
Kahler potential in terms of a meson M has been determined in Ref. [JT], 

K = ^logTrVMMt . (5.18) 
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If we relax the vacuum moduli to be equal ( 15.41) . thus having the possibility of distinct {uj}'s in 
Eq. (15. ip . it will prove convenient to work directly with the mesons of SO(2M) 

(e u (z — a) ±ia \ 
, 5.19 
±ia e~ u (z + a)J 

with a,u G M. The meson VEV will be diag(uf , uf) = diag(e u , e~ u ). Using the Kahler potential 
(15.181) we readily obtain the energy density 

E = £dd log ( \z - a tanh(w) | 2 + %— ] . (5.20) 

V cosh (u) J 

Furthermore, we can construct a typical example of fractional vortices, in a U(l) x SO(2N) 
model in the lump limit as follows 

E = 2idd log |2 ~ a » tanh (^)| 2 + cq 1 J , (5-21) 



,i=l 



with ufj-i = rriie Ui ,V2i = rriie~ Ui . For each 50(2) subgroup, we have in this construction 
a possibility for amplification and a^Wj which serve as a position- and an effective size- 
parameters. One can observe that m, controls the relative weight of the energy distributed to 
the i-th fractional vortex. 



6 Conclusion 

In this paper we have given a simple account of the fractional vortices which have minimally 
quantized magnetic flux (winding) but with non-trivial substructures in the energy distribution 
in the transverse plane. They could often appear in various generalizations of the Abelian Higgs 
model. The common characteristic features these models share are a non-trivial vacuum degen- 
eracy and the BPS saturated nature of the vortex solutions. We have generalized the moduli 
matrix formalism [201 |2~T1 [22] to clarify all possible moduli parameters of the minimally quantized 
fractional vortices. 

The vacuum moduli M. in these models turns out, in general, to be a singular manifold, i.e., 
a manifold with singularities. Vortex solutions approach the vacuum configuration far from the 
center, and trace various closed gauge orbits F. 

We have classified fractional vortices into two types; the first type appears when Ai has a 
Z n singularity where the gauge symmetry is not restored while the second type occurs when 
M. has a 2-cycle with a deformed geometry. The existence of a singularity is not essential for 
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the fractional vortices of the second type. Indeed, we have observed that smooth fractional 
lump solutions become singular as the smooth manifold Ai is deformed into a singular manifold 
(e.g. when some FI parameters are turned off). Even when Ai has such singularities, we have 
found smooth fractional vortex solutions. The vortices share the same properties as those of the 
corresponding lumps wrapping on M. smoothened. 

An interesting aspect of our analysis, especially relevant to the systems with the first type 
of fractional vortices and lumps, is the fact that the latter often represent a generalized fiber 
bundle over the singular manifold M. . At a singularity (or on a singular surface) the fiber space 
F undergoes a discontinuous change either in its dimension or in its nature, or both. The vortex 
moduli also make a jump at such a point (points). These observations seem to point towards 
interesting physical applications as well as some novel kind of mathematical structures. 



Note added 

After completing this paper, a new paper by D. Tong and B. Collie which discusses precisely the 
second type of our fractional vortices - lumps, though in a different context, was posted on the 
ArXiv El. 
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A General BPS vortex equations 



A.l General arguments 

Let denote the gauge fields for G = U{l) n x G' gauge group, where 7 = — (n — 1), • ■ ■ , — 1, 
are for the Abelian gauge fields ?7(l) n , the rest referring to the arbitrary simple group G' . We 
write the complex scalar (Higgs) fields as 

H a i a= (r 1 ,A 1 ) ,(r 2 ,A 2 ) 

where T{ and A\ stand for the color and flavor indices of the i-th matter in the representation R{. 
A gauge transformation (G) induces 

5H a = i (AH) a = i k \ T % H a = iJ2 Al ( T % H^ + iJ^ Al ( Tl )p H ^ , 

I I<0 I>0 

Here T 1 for I < is the charge matrix for U(l)i, i.e, Q{ are the U(l)r charges of the matter 
fields present 

{T% = Q 1 ^^ , for 7<0. 

T 1 (I > 1) denotes the generators of the gauge group G' in the (in general, reducible) represen- 
tation 

T 1 = tf ] ® 1^(1) © tf ® 1^(2) © ■ ■ ■ . 

Accordingly, the covariant derivative is given by 

(V,HY = [(d, - iW») Hf = [6p, - iWl (T%) ; 
in our models the energy (tension) has the following form and is semi-positive definite: 

J l i 9i a j 

where f j = 0, ( v 7 > 0) and £/ for I < are assumed such that a supersymmetric vacuum exists 
and the theory is in Higgs phase. 

The vacuum manifold is the zero locus of the superpotential 

M v = {H | D 1 = WT 1 H -6 = 0, V/} (A.l) 
and the vacuum moduli is 

M v = M/G ~ H//G € . (A.2) 
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These quotients may be ill-defined where the vacuum moduli space is singular or when some 
gauge symmetry is restored. 

When we choose a point in the vacuum manifold where some of U(l) gauge symmetries are 
spontaneously broken as G = U(l) n x G' — »■ U(i) n ~ r x g' (g' E G'), the topologically stable 
vortices appear with support from the non-trivial first homotopy group 

r 

8=1 

This is related to n 2 (Ai) by the homotopy sequence (I3.8p . We now introduce the complex 
coordinates to write down the BPS equations 

1 _1 
-(d x -id y ), d=-\ 



z = x + iy, z = x-iy, d = -(d x - id y ) , d = -(d x + id y ) , 



V = i(X>! + iV 2 ) , V = 1 -{V X - iV 2 ) , W = ^{Wl - iWl)T l , w = i(w/ + %W l 2 )T l . 
The Bogomol'nyi completion reads 

T = [d 2 x T g^T 1 !! - 0)] 2 + \( VlH ± ^#)T =F CiFL 

This shows that the minimum of the tension is given by either the "negative chirality" (or left- 
winding) solutions satisfying the Bogomol'nyi equations 

F( 2 = gj^T'H - fr) , VH — , T min = J d 2 x F* 2 > ; (A.3) 
or by the "positive chirality" (right-winding) solutions such that 

F[ 2 = -g&rfT'H - £j) , VH = , T min = 6 ^ rf 2 x F( 2 > ; 
Therefore, the tension is expressed by 

T min = 2 <i v i > Ul = T ^T I dx " , (A.4) 



where — is for the negative chirality and + is for the positive chirality. The rational number vj 
(for / < 0) stands for how many times the corresponding solution winds U(l)i ~ S 1 when we go 
around once the boundary circle S 1 on the z-plane. 

From now on, we concentrate on the BPS states with the negative chirality which are the 
solutions of Eq. flA.3j) . The matter part of the above vortex equation can be solved by the moduli 
matrix H Q (z) which is a color-flavor mixed matrix whose elements are holomorphic (polynomials) 
in z, for the left-handed vortex solution Eq. (TA~3!) [20112X1122] 

H = S~ l H {z) , W = -iS^dS , (A.5) 
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where S is an element of the complex extension of the gauge group 

S{z, z) = e< E ' kI{z ~ z)T ' e G c , A J (z, z)eC. (A.6) 

For a given Hq, S is determined by the first equation in Eq. (1A.3|) . We assume existence and 
uniqueness of the solutions. Thus all the complex constants (coefficients of the polynomial 
functions) appearing in Hq(z) represent the vortex moduli parameters. Finally, S and Hq(z) are 
defined up to a V(z) transformation which does not change (1A.5|) and keeps S inside G c 

H (z)^V(z)H (z) , S^V(z)S , 

where V(z) is a holomorphic matrix belonging to the complexified gauge group G c . When the 
vortex moduli space is seen as a complex manifold (whose local coordinates are the moduli 
parameters appearing in Hq(z)), the V(z) transformations act as the transition functions in 
two overlapping patches. Thus the (vortex) moduli space of the 1/2 BPS vortices is formally 
expressed by 

M vov = {H | H ~ VH , BH = , dV = 0, V e G c } . 

Generally speaking, one should choose a boundary condition when one wishes to solve some 
partial differential equations. Our strategy of solving the BPS equations (IA.3I) is somehow 
upside-down to such an ordinary way. In fact, we have not fixed any boundary conditions yet. 
Of course, we are talking about the topological solitons which are characterized by the boundary 
conditions. So the remaining task is to figure out the condition for the moduli matrix which yields 
solutions satisfying the correct boundary condition. Note that the condition for Ho(z) depends 
on the boundary condition which is nothing but VEV (H), a point on the vacuum manifold. 
We have to be careful to specify the moduli matrix especially when we choose a singular point 
as the boundary condition. Furthermore, the configuration (energy distribution) may change as 
varying the VEV even if the moduli matrix is fixed. 

The concrete conditions for the moduli matrix H (z) depend on details of the models, such as 
gauge groups, representations of the matter fields and the U (1) charges. The case of G = U(l)xG' 
with G' being an arbitrary simple group have been studied in Ref. [38]. In what follow, we will 
explain two typical cases, i) G = U(l) gauge theory with matter fields whose ?7(l)-charges are 
distinct, ii) there are more than one Abelian group such as G = U(l)i x f/(l) 2 . These models, 
and especially, their moduli spaces, have not been studied so far. 
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A.2 G = U(l)i x 17(1)2 



Let us consider a, G = U(l)i x {7(1) 2 gauge theory with three Higgs fields if = (A,B,C) with 
the following {7(l)-charges: Qi = (m, 1, 1) under the first U(l)i and Q 2 = (0, 1, —1) under Z7(l) 2 . 
The vacuum manifold and the vacuum moduli are 

M = {A,B,C\m\A\ 2 + \B\ 2 + \C\ 2 = ^, \B\ 2 -\C\ 2 = ^ 2 }^S 3 , (A.7) 
M = M/({7(l)i x C/(l) 2 ) ~ CP x /Z m . (A.8) 

Here £1,2 are the Fl-terms of {7(1)1,2- We choose them in the region — £1 < £ 2 < £1 (£1 > 0) to get 
the system in a Higgs vacuum. In the main text of this paper we set £ 2 = but here we consider 
a more general situation. At generic point of M, both £7(l)i and {7(1) 2 are broken, hence we 
have topologically stable vortex solutions with topological characters uj 

Vl = ^- J dx 2 dd\og\ Sl \ 2 , W^-idlogsr, (7 = 1,2), (A.9) 
H = (A, B, C) = (s^ m A (z), s^s^Boiz), s^s 2 C (z)) , (A.10) 



The solutions are determined by C7 (l)i x U(l) 2 winding numbers (z/i,z/ 2 ). Eq. ( 1A.9I) determines 
the asymptotic behavior of sj 

\si\ 2 -> |z| 2 " J as |z| -> 00 . (A.ll) 
Then the Higgs fields asymptotically behaves as 

{A ,B Q ,C )~(\z\ m ^A vev , \z\^B vcv , \z\^C vcv ) , as |z| - 00 , (A.12) 
(A vev , 7? vcv , C vcv ) being a point in the vacuum manifold. 

A. 2.1 m = l 

Let us first consider a generic point, i.e., (A vev , B vev , C vcv ) 7^ (0,0,0). Holomorphy forces us to 
choose V\ G Z + , z/i + z/ 2 G Z + and i>i — i> 2 G Z + . Let us rewrite 

1/1 = fci G Z + , ui + ^2 = A; 2 G Z+ , 2fci > A; 2 > . (A.13) 

The tension can be expressed as 

L = + V £ 2 = fa - ^k, + &k 2 > > . (A.14) 

Z7T Z 

Thus the moduli matrix for (k\, k 2 ) configuration is 

A (z) = A vcv z k > + • • • , B (z) = B vcv z k * + • • • , C (z) = C vev z 2k ^ + ■■■ . (A.15) 
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The number of complex moduli parameters is dimcAi(k l ,k 2 ) = The minimum configuration 
depends on £ 2 - When < £ 2 < £1, the minimum configuration is T^ q) = 27r(£i — £ 2 )- When 
— £1 < £2 < 0, T(i t 2) = 27r(^i + £ 2 ) is minimum. The next lightest configuration has always 
T^i) = 27r£i. If £2 = 0, all these three are degenerate, namely T(i j0 ) = = ^(1,2) ■ 

Let us next choose the case of A vcv = with B vev , C vev 7^ 0. In this case, the condition 
changes as v\ + v 2 = k\ G Z + and v\ — v 2 = k 2 G Z+. In this case, the ?7(l)-charges ^2 are half 
quantized and the energy is 

_ = -^-A* + _^_Jfc 2 > . (A.16) 

The moduli matrix for (&i, fo) configuration should be chosen as 

A (z) = az a + ■ ■ ■ , B (z) = B vcv z k ^ + • • • , C (z) = C vcv z k * + ■ ■ • , (A.17) 

where a is an arbitrary positive integer less than V\ = (A;i+A; 2 )/2. The number of complex moduli 
parameters is dimc-M^^) = k\ + k 2 + a + 1. The minimum configuration is Tn i0 ) = 7r(£i + £2) 
for — £1 < £ 2 < while T(o,i) = ?r(£i — £2) for < £2 < Ci- ^(1,1) = 2?r£i is the third lightest 
configuration. Note that these tensions are one half of that of the vortices for A vcv 7^ 0. If £2 = 0, 
we again observe degeneracy TVq) = T(o,i) < ^(1,1) • 

The vacuum _B vev = C vcv = is possible when £ 2 = 0. We shall not consider this case since 
the second U(l) factor is now in a Coulomb phase. 

A.2.2 m = 2 

Let us repeat the analysis for case m = 2. First we consider (A vcv , B vcv , C vcv ) 7^ (0, 0, 0). We find 
the condition 2v\ G Z+, V\ + i/ 2 G Z + and 1/1 - ^ G Z + . Rewrite 

^ + v 2 = k\ G Z + , 1/1 — ^2 = k 2 G Z + . (A. 18) 

Note that 2v\ = k\ + k 2 G Z + is automatically satisfied. The energy is given by 

^ = ^1^— + fc-g— = + ~^ k * > • (A.19) 

Since the U(l) charges are half quantized, this tension formula is slightly different from Eq. (1A.14I) 
with m — 1. The moduli matrix for (ki, k 2 ) configuration is 

A (z) = A vcv z k ^ + • • • , B (z) = B vcv z k ^ + • • • , C Q (z) = C vcv z k * + ■■■ . (A.20) 

Dimension of the moduli space is dimc-Mo^^) = 2k\ + 2k 2 . The minimal energy configuration is 
^(1,0) = ^(^i — £2) for £2 > or T(i 1) = 7r(£i + £2) for £2 < 0. If £2 = 0, these two are degenerate. 
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Configurations with A vcv = 0, or B vcv = (or C vcv = 0) as the boundary condition are not 
special but belong to the above category. B vcv = C vcv = is only possible when £ 2 = 0. Again, 
in this case the (1)2 gauge symmetry is restored (Coulomb phase). 

A.2.3 m>3 

As we have seen, the difference coming from the choice of m appears only when (A vev , B vev , C vev ) 7^ 
(0,0,0): 

2 Ah 

mv\ = ki , V\ + V2 = &2 j =>• i>i — b>2 = 1v\ — &2 = *2 £ Z + . (A. 21) 

m 

When m is odd, we should choose k\ = m odd /^ with k[ G Z+. Then k 2 can be taken in the region 
2k[ > k2 > 0. The corresponding tension is 



^odd _ £1 ~ £2 
2tt m odd 

The moduli matrix is 



ki + 6*2 = (6-6)*'i + £2*2 ■ (A.22) 



Aq(z) = A vcv ^ fc 'i + • ■ ■ , fl (z) = £ vcv ^ 2 + ■ ■ • , C (z) = a c v^ 2fci - fc2 + • • ■ • (A.23) 
The dimension of the moduli space is dim c A^°^ i d fc2 ^ ) = (m odd + 2)k[. 

When m is even, we must choose k\ = m ^ en k[ with k[ G Z+. Thus the C/(l) charges are half 
quantized. Holomorphy forces k[ > k 2 > 0. The tension is given by 

= ^^*i + 6*2 = ^^*i + 6*2 • (A.24) 

Z7T Ul even Z 

The moduli matrix is given by 

A (z) = A vev z^^ + • • • , B (z) = B vev z h > + ■■■ , C (z) = C vcv z k '^ + ■■■ . (A.25) 
The dimension of the moduli space is dimcA^^ n fc2 ^ = ( m< ^ cn + l) k[. 

A.3 G = U{1) x SU(N) 

Let us consider G = U(l) x G' and G" = SU(N) with Higgs fields if = (A, B) where A is 
a singlet scalar field and B is a collection of iV fields in the fundamental representation N_ of 
SU(N), written as an N x N matrix. We assign the U(l) charges (m, 1) to the fields (A, B), 
respectively. The vacuum manifold is given by the D-flatness conditions and the vacuum moduli 
is obtained by dividing it by G as 

M= {(i,B) I m\A\ 2 l N + NBB^ = £1 N } , A* = M/G ~ CP 1 /Z m .. (A.26) 
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The vacuum condition is simplified by using B = B vcv 1n as 

m\A vcv \ 2 + N\B VCV \ 2 =£ . (A.27) 

The moduli matrix formalism is summarized as 

(A,B) = s' 1 (A (z), S^B (z)) , Wu(i) = -tdlogs , W SU (n) = -iS^dS , (A.28) 

with s G C* and S G SL(N, C). U(l) winding number is given by 

v = ^- J dx 2 dd log |s| 2 , \s\ 2 ~\z\ 2u as \z\ -»• oo . (A.29) 

When we have some non-Abelian gauge group G', we can consider the holomorphic G' invariant 

I = det B = s~ N det B (z) -> I*! - *" det B (z) = B^ cv as |z| -»■ oo . (A.30) 

Since all the elements in Bo(.z) is holomorphic in z, det Bq(z) is also a polynomial function of z. 
Combining this and the behavior of A 

A = s _1 A,0) -> |^|- m!/ A (2) = A vev , as |*| -»■ oo , (A.31) 

we can consistently determine v. 

When we choose (A vcv , B vcv ) ^ (0,0), the moduli matrix asymptotically behave as 

A (z) -> \z\ mv A YeY , detB (^) -> laf^B^ , as |z| -+ oo . (A.32) 

Thus we should choose the U(l) winding number v in such a way that 

jV N' 

mu = keZ + , Nu = —k = —keZ + , (A.33) 

m m 

where we have assumed N = n N' an d m = rio^n' with n G Z + (no = g.c.d(A, m)). Thus fc must 
be = m'k' with fc' G Z+, so that the t/(l) winding number is fractionally quantized v = ^ 
if no > 2. The corresponding tension is T = 27r£— . The moduli matrices for the configuration 
with k l may be chosen as 

A (z) = A vev2 ;"o fc + ■ ■ ■ , det flo(z) = ^ v ^^ + ■ • • • (A.34) 
Note that we have specified only det ~Bq(z), so there are many possibilities for B (z) per se, it is 

N+m 

nothing but the orientational moduli. The moduli space for k — 1 is M.y=\ ~ C ™o x CP . 

When we choose 74 vev = and B vcv ^ as the boundary condition, we should satisfy 

A -> |z|- m M (z) = , / -> 1^1"^ det B (z) = as |z| -»■ oo . (A.35) 

Thus we can choose Nu = kE Z + , such that the U(l) winding is 1/N quantized and correspond- 
ing tension is T = 2ir£j7. The moduli matrices are 

A {z) = az a + ■ ■ ■ , det B = B^ v z k + • ■ ■ , (A.36) 

where a is a semi-positive definite integer a < mu = mjfr. 
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